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Abstract. We present a calculation of the elliptic flow and azimuthal dependence of the correlation radii 
in the ellipsoidally symmetric generalization of the Buda-Lund model. The elliptic flow is shown to depend 
only on the flow anisotropy while in case of correlation radii both flow and space anisotropy play an im- 
portant role in determining their azimuthal oscillation. We also outline a simple procedure for determining 
the parameters of the model from data. 
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! 1 Introduction 

Studies of the freeze-out configuration in non-central ultra- 
. relativistic nuclear collisions have yielded valuable infor- 
' mation on the dynamics of the collision. The azimuthal 
, asymmetry of hadronic momentum spectra — widely known 
■ under the somewhat vague term "elliptic flow" — is a cor- 
I nerstone result for conclusions about thcrmalization and 
. near perfect fluid dynamics [1]. The measurement of az- 
imuthal dependence of correlation radii in non-central col- 
lisions indicates a fireball which breaks up while still be- 
1 ing elongated out-of-reaction-plane [2] . This is the original 
orientation of the system and the observation thus pro- 
vides an upper bound on the total lifespan of the created 
fireball and suggests an early and sudden freeze-out in a 
blast-wave picture, influentially suggested in ref. |3j. 

These observables can result from an interplay of spa- 
tial asymmetry of the fireball and the asymmetry of the 
flow velocity field of an exploding fireball. The question 
is if and how these asymmetries can be both determined 
from data. Such a question can be studied in the frame- 
work of a model. Here we present a study that utilizes the 
Buda-Lund hydro model. 

For the present study we use the ellipsoidal Buda-Lund 
parametrization of the freeze-out of the fireball [4]. Note 
that the parameterization corresponds to a solution of cer- 
tain class of hydrodynamic models [5l6,7,8j. 

We shall show that in the Buda-Lund model elliptic 
flow results solely from the flow asymmetry, and that the 
entanglement of spatial and flow anisotropy in determin- 
ing the azimuthal oscillation of the correlation radii is 



rather strong. This can be contrasted to other derivations 
that explore the conditions for quark coalescence [9], or 
the numerical solutions of boost-invariant, azimuthally de- 
pendent hydrodynamics with Gaussian initial conditions 
jlOj . Azimuthally asymmetric HBT radii were also con- 
sidered recently in cascade models, e.g. in the fast Monte- 
Carlo model of ref. [Tl], or, in the Hadronic Resonance 
Cascade [H]. A simple procedure to disentangle spatial 
and momentum space anisotropy from elliptic flow and 
azimuthally sensitive HBT measurements was proposed 
for generalized blast-wave type of models in reL 13J. Now 
we investigate how to disentangle these effects in the ellip- 
soidally symmetric version of the Buda-Lund hydro model. 
Let us note that we cannot attempt to evaluate or fully 
cite all the relevant works in this technical paper, but we 
kindly refer the interested readers to the review of reviews 
in ref. '14], and also to ref. [15], a recent review that em- 
phasizes the femtoscopy aspects of particle interferometry. 



The structure of our paper is as follows: In Section [5] 
we shall shortly introduce the basic features of the el- 
lipsoidally symmetric Buda-Lund model. Section [3] deals 
with the calculation of elliptic flow and Section 2] with de- 
termining the correlation radii. Finally, in Section [5] we 
present a simple strategy how various parameters of an 
azimuthally non-symmetric fireball can be obtained from 
data. We conclude in Section [6l 
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2 Buda-Lund model: basic features 

Perfect fluid hydrodynamics is based on local conservation 
of entropy a and four- momentum tensor T^"^ , 



^^{au^^) = 0, 

a^T^'' = 0, 



(1) 

(2) 



where it^ stands for the four-velocity of the matter, which 
is normalized to unity as u^^u^ — 1. (Entropy conservation 
indicates that there are no dissipative terms like bulk or 
shear viscosity, and heat conduction, as any of these effects 
would lead to entropy production). For perfect fluids, the 
four-momentum tensor is diagonal in the local rest frame. 



T^'' = (e + p)u^'u'' - pg 



fJ.V 



(3) 



Here e stands for the local energy density and p for the 
pressure. The local conservation of energy, momentum and 
entropy yield five constraints on six unknowns (the three 
independent components of the four- velocity field and the 
three thermodynamical variables, the energy density, the 
entropy density and the pressure). These equations are 
closed by the equation(s) of state, which gives the rela- 
tionship between e, p (and a), sometimes as a function of 
yet another variables, like the temperature T. 

We focus here on the analytic approach in exploring 
the consequences of the presence of such perfect fluids in 
high energy heavy ion experiments in Au-|-Au collisions 
at RHIC. Such exact analytic solutions were published re- 
cently in refs. [5l[6ll7l[8jll6j. In some of the recently found 
exact analytic solutions of relativistic hydrodynamics, the 
equation of state has a relatively simple form, typically 
e — B — k{p + B) is assumed, where i? is a constant, called 
the bag constant which has a non-zero value in a quark- 
gluon plasma and has a vanishing value in the hadron gas 
phase. The parameter k is either a constant like in ref. [16] 
or more generally, analytic solutions can be obtained even 
if K = k(2^) is an arbitrary, temperature dependent func- 
tion, see for example the non-relativistic exact solutions 
in ref. [5]. 

A tool, that includes these above listed exact, dynam- 
ical hydro solutions, as well as some other cases, and in- 
terpolates among them, is the Buda-Lund hydro model of 
refs. [TTlfTS]. This hydro model is successful in describ- 
ing experimental data on single particle spectra and two- 
particle correlations [4lll9|. 

The Buda-Lund hydro model [17] successfully describes 
BRAHMS, PHENIX, PHOBOS and STAR data on iden- 
tified single particle spectra and the transverse mass de- 
pendent Bose-Einstein or HBT radii as well as the pseu- 
dorapidity distribution of charged particles in central Au 
+ Au collisions both at y^SNN = 130 GeV [19] and at 
Y^snn = 200 GeV H] and in p-l-p colhsions at -^/i = 200 
GeV |2D], as well as data from Pb-^Pb colhsions at CERN 
SPS [H] and /i -hp reactions at CERN SPS |22l|23]. The 
model is defined with the help of its emission function; to 
take into account the effects of long-lived resonances, it 
utilizes the core-halo model 1241. 



In the Buda-Lund model, the emission function is given 
by that of a hydrodynamically expanding fireball, sur- 
rounded by a halo of long lived resonances. In this paper, 
we assume the validity of the core-halo model, i.e. that the 
long-lived resonances create an unresolved, narrow peak in 
the two-particle HBT correlation function, which results 
in an effective intercept parameter A. The integrals of the 
emission function of the hydrodynamically expanding core 
are evaluated |22| in a saddle-point approximation. When 
an improved saddle-point calculation is performed, effec- 
tive binary sources and related non-Gaussian behavior is 
seen in this model, in particular in the longitudinal direc- 
tion [22j . However, in a Gaussian approximation, using a 
single saddle point in the integration, the effective Buda- 
Lund emission function looks like: 



S{x, k)d'^x 



g pf'u^{xf)H{Tf) 



(4) 



(27r)3 B{xf,p) + Sq 
X exp [~R'^'^{x - XfY{x - XfY] dt^x, 

where g is the degeneracy factor [g — 1 for identified pseu- 
doscalar mesons, g = 2 for identified spin=l/2 baryons), 
B{x,p) is the (inverse) Boltzmann phase-space distribu- 
tion, the term Sq is determined by quantum statistics, Sq = 
0, —1, and -|-1 for Boltzmann, Bose-Einstein and Fermi- 
Dirac distributions, respectively, and the space-time four- 
vector of the point of maximal emissivity (or freeze-out) 
Xf = Xf{p) depends on the momentum of the emitted par- 
ticles. The time dependence of the emission is described 
by H{t), and 



df,d^{-\TL{So)){xf,p), 



(5) 



is the 'width' of the emission function, introducing 5*0, as 
the 'main' part of the emission function: 



So{x,p) 



H{t) 



B{x,p) + Sq 



(6) 



For a relativistic, hydrodynamically expanding system, 
the (inverse) Boltzmann phase-space distribution is 



B(x,p) — exp 



p ■ u{x) fJ,{x) 
T{x) ~T(x) 



(7) 



We will utilize some ansatz for the shape of the flow four- 
velocity, Ui,{x), chemical potential, /i(a;), and temperature, 
T{x) distributions. Their form is determined with the help 
of exact solutions of hydrodynamics, both in the relativis- 
tic [6,16,8J and in the non-relativistic cases ^25,26.27.,28]. 
with the conditions that these distributions are character- 
ized by mean values and variances, and that they lead to 
(simple) analytic formulas when evaluating particle spec- 
tra and two-particle correlations. 

The solution is scale-invariant with the following scal- 
ing variable: 



X2 ^ y2 



Z2' 



The velocity field is 



(8) 



(9) 
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with 7 being the appropriate factor to ensure u ■ u ~ 1. 

For the fugacity distribution we assume a shape, that 
leads to Gaussian profile in the non-relativistic limit, 



T{x) 



To 



(10) 



corresponding to the solution discussed in refs. [25ll26l l5|. 
We assume that the temperature may depend on the posi- 
tion as well as on proper-time. We characterize the inverse 
temperature distribution similarly to the shape used in the 
axially symmetric model of ref. [171129] . and discussed in 
the exact hydro solutions of refs [ISIIS], 



1 



Tix) 



1 



Tq-Ts 



To - Te (r - Tof 



(11) 

where Tq and are the temperatures of the center, and 
the surface at the mean freeze-out time tq, while Tg cor- 
responds to the temperature of the center after most of 
the particle emission is over (cooling due to evaporation 
and expansion). Sudden emission corresponds to Tg = Tg, 
and the At — )■ limit. It is convenient to introduce the 
following quantities 



2 ^0-^^ /^T\ 
,2_ To-Te _ /AT\ 



(12) 
(13) 



The time dependence of the emission, described by 
H{t) can be approximated with a Gaussian representa- 
tion of the Dirac-delta distribution around the freeze-out 
proper-time r/. 



H{t) 



1 



(27rZ\T2)i/2 



exp 



2Z\t2 



(14) 



with At being the proper-time duration of the particle 
production. 

From now on our investigations will be performed at 
midrapidity. With the above shorthands, at niidrapidity 
the Boltzmann factor at the freeze-out will have the fol- 
lowing form: 



B{xf,p) = exp 



pI . pI Pt , fio 



2mtTx 2mtTy 2mtTo Tq Tq 



(15) 

where mt = -f and the direction dependent slope 
parameters are 



=To + mtX 



2 To 



To + mto^ ' 



Ty =To + mt Y 



>2 ^0 



To -I- mto? 



(16) 
(17) 



3 Elliptic flow 

The result for the elliptic flow, that comes directly from 
a perfect hydro solution is the following simple scaling 



law 



V2 



hjw) 



(18) 



where /„(z) stands for the modified Bessel function of 
the second kind, In{z) = (I/ti") /J^ exp(z cos(0)) cos(n6')(i6'. 
The scaling variable 



/ 1 1 



4mt \Ty Tx 
This can also be written as 



w = E 



K 



T 



cff 



(19) 



(20) 



where Ek is a relativistic generalization of the transverse 
kinetic energy, defined as 



E 



K 



2mf 



(21) 



and we have introduced Teff, the effective slope of the az- 
imuthally averaged single particle pt spectra as the har- 
monic mean of the slope parameters in the in-plane and 
in the out-of-plane transverse directions, 



Tnff — 2 



1^1 



and momentum space eccentricity parameter 



Tx ~h Ty 



(22) 



(23) 



If there is no transverse temperature gradient (a = 0), 
then this simplifies as 



x2_r2 



(24) 



(see refs. [511151150] for details). 

Thus the Buda-Lund hydro model predicted [51IT8] a 
universal scaling of the elliptic flow: every V2 measurement 
is predicted to fall on the same scaling curve Ii/Iq when 
plotted against the scaling variable w. This means, that V2 
depends on any physical parameter (transverse or longi- 
tudinal momentum, mass, center of mass energy, collision 
centrality, type of the colliding nucleus etc.) only through 
w. 

We note, that for a leading order approximation, Ek ~ 
mt — m, which also explains recent development on scaling 
properties of V2 by the PHENIX experiment at midrapid- 
ity [3I1|32J. 

In general, the azimuthal anisotropy of hadronic mo- 
mentum spectrum is generated by the azimuthal variation 
of both the absolute magnitude of the expansion velocity 
and its orientation. The same elliptic flow parameter V2 
can be obtained by different combinations of these two 
and it has been shown in the framework of blast-wave 
model that by modifying the latter one can completely 
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change (even invert) the dependence of V2 on the for- 
mer [33]. TechnicaUy, this shows up as a dependence of 
the eUiptic flow on both the spatial anisotropy and the 
flow anisotropy. 



We see from eqs. ([T8 | - ([24l) that in Buda-Lund model 
the velocity fleld is oriented in such a way that the spatial 
anisotropy of the hadronic freeze-out stage plays no role 
in determining V2- eq. (|18p depends only on the transverse 
expansion rates X and Y , but it does not depend on the 
actual source sizes X and Y. Only difference of the slope 
parameters and Ty is important for the elliptic flow. 
Thus the elliptic flow vanishes in this model, ii X = Y, 
indicating that if the expansion rates are similar in the 
in-plane and out-of-plane transverse directions, V2 = 0. 
In this sense, we may write that V2 in the Buda-Lund 
model is driven by the difference of the in-plane and out- 
of-plane transverse expansion rates of the fluid, during the 
hadronic final state. 

In order to compare more easily the Buda-Lund model 
results for the elliptic flow with the azimuthally sensitive 
extension of the blast- wave model ;2,,33, , we introduce po 
and p2 so that 



X = Po{l+p2) 

Y = poil-p2) 



(25a) 
(25b) 






-0.6 -0.4 -0.2 0.2 0.4 0.6 
P2 



1 

0.8 
0.6 
0.4 
0.2 





-0.6 -0.4 -0.2 0.2 0.4 0.6 
P2 



therefore 



Po 



P2 = 



2 

X -Y 
X + Y 



(26a) 
(26b) 



The elliptic flow in the Buda-Lund model also depends 
on the transverse temperature gradient a of eq. (jl2[) . It 
actually moderates the difference between and Ty and 
increases the elliptic flow at given p2. The dependence of 

on flow anisotropy p2 and temperature gradient a is 
shown in Fig. [T] 



4 Azimuthally dependent correlation radii 



Fig. 1. Lines of constant V2 of pions as function of p2 and a. 
The increment between neighboring contours is 0.05, the thick 
curve shows V2 — 0. The sign of V2 is the same as the sign 
of p2. Transverse momentum is 0.5 GeV/c (top) and 1 GeV/c 
(bottom). Other parameters in the calculation were motivated 
by data (see Section[5l): To = 163 MeV, po = 0.61. 



with A being the intercept parameter, q = pi ^ p2, and 



^ — X 



1 mt (P 



Hi 



Tn 



;, '"-(°^+^') v'. 



(28a) 
(28b) 

(28c) 

(28d) 



The two-particle Bose-Einstein correlation function from 
the Buda-Lund source function (defined in eq. (jl])) is most 
naturally calculated in the frame given by the main axes 
of the ellipsoid. The result is 



(27) 



Note that parameter d makes the effective life-time and 
the difference between the out and the side radius compo- 
nents momentum dependent. 

From the mass-shell constraint of p\ 
finds that 



m one 



(29) 
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where /3 = K / K'^ = (pi — P2) / {pi — P2) ■ Thus we can 
rewrite eq. ([27|) to 



C7((z) = l + Acxp - ^llilj] (30) 

with the modified radii of 

Rl = RX^ + 0iATl, (31a) 

Rl^RXy^HlArl, (31b) 

rI = RX,+151AtI, (31c) 

i?L = /3x/3.^t2, (31d) 

Rl, = PyP.Arl (31e) 

Rly = P.PyArl (31f) 

The above formulas are given in the frame, specified by the 
main axes of the expanding, ellipsoidally symmetric fire- 
ball. In general, this frame is rotated [5] with respect to 
the standard out-side- long system [34J , which is given by 
the beam direction (longitudinal) and the transverse com- 
ponent of the pair momentum (outward). The correlation 
function in this out-side-long system is parameterized as 

C{q) = I + Xc-'JoRl-liRl-lfR^-^i']oqsRl,+qoqiBli+q,q,R^^l) ^ 

(32) 

We denote by i? the tilt angle between the longitudinal 
direction and the z-axis and ip stands for the angle be- 
tween the outward direction and the plane given by x and 
z axes — the reaction plane. In the new frame we obtain 
azimuthally sensitive HBT radii that are given as 



Rl = R'^ cos^ ^ + Rl sin^ ip + R'^y sin(2(^) , (33a) 
Rl = R'^ sin^ ^ + Rl cos^ ip - R% sm.{2ip) , (33b) 
Rf = Rl sin^ 1^ + Rl cos^ ^ + Rl, sin(2i9) , (33c) 
2Rl, = -Rl sin(2(^) + i?^ sin(2^) + 2R'^y cos(2(^), (33d) 
2Rli = {Rl sin(2i9) + Rj sin(2i9) - 27?^^ cos(2i9)) sin 9? 

+ {2RlySmd + 2Rl,cosd) cosip, (33e) 
2Rl = [Rl sin(2z9) - R^ sin(2i?) -f- 2R^, cos(2z9)) cos(^ 
+ {2Rly sin 1} + 2Rl, cos i?) sin ip, (33f ) 

where we have introduced 

R'^ = Rl cos^ i) + Rl sin^ d - Rl, sin(2z9) (34) 
R'xy = Rly cos d - RI, siu 1? (35) 



If the fireball is not tilted in the reaction plane (i? = 0), 
which might appear at very high collision energies, then 
in the longitudinally co- moving frame (LCMS, /3; = 0) the 



formulas further simplify as 

Rl = i?L cos^ + Rl y sin^ ^ + PIAtI (36a) 

n2 I C>2 p2 _ 7->2 

^*.x ' ^*.y , o2 A 2 ^*.y ^*.x /r, \ 

= — ^ '- + PoAt: ^ cos(2v3) 

Rl = Rl ^ sin^ ip + Rl y cos^ ip (36b) 

f>2 I E>2 d2 _ p2 

= — ^ ^ + cos(2^) , 

i?L = ^^^-^Bin(2^), (36c) 

i?o;=0, (36e) 
i?'i-0. (36f) 



It is convenient to decompose the correlation radii into 
Fourier series in azimuthal angle [351136] . Here, it is partic- 
ularly trivial since there are only zeroth and second order 
terms 

Rl = Rl^ + Rl2Cos{2^) (37a) 

Rl^Rl^ + Rl2Cos{2^) (37b) 

Rl = Rh (37c) 

C -Co + C2sin(2^)- (37d) 

In this Gaussian single saddle point version of the 
Buda-Lund model there is no implicit azimuthal depen- 
dence of the homogeneity regions. Thus the complete de- 
pendence on ip is written in eqs. ([36]). In such a case we 
observe rather known features: 

- The oscillation amplitudes 2I' l-^s 2I' ^^cl \Rlg 2 1 
equal. 

- The effective duration of particle emission can be ob- 
tained from the difference of non-oscillating parts as 

PIAtI - i?2 - Rl^ ■ (38) 

- We can extract the effective source sizes Rl^^ and Rl y 
from data as 

RI,x = Rlfi - R^s,2i (39a) 
Ri,y = RI,q + RI,2 ■ (39b) 

If the first feature is not observed, or, if data lead to a neg- 
ative r.h.s. of eq. ([55]) . this would indicate that the above 
presented Gaussian version of the Buda-Lund model is 
not applicable for the description of that given data. It is 
however well known, already in the first formulation of the 
Buda-Lund model, ref. [17], the r.h.s. of eq. ([38]) can be 
negative. The exact hydro solutions and symmetry argu- 
ments require that Rsfi = Ro.o in the vanishing transverse 
momentum limit only, otherwise at finite transverse mo- 
mentum, Rs Q > Rofi is allowed in general. Such a behav- 
ior happens in exact, ring of fire type of hydrodynamical 
solutions, which were found in fireball hydrodynamics in 
ref. [55]. Such rings of fire type of exact hydro solutions 
appear if the radial flow is moderate and the temperature 
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distribution is very inhomogeneous, in our present nota- 
tion if c? ^ m,p\l (2Tq), while nearly Gaussian fireballs 
appear, if the temperature is nearly homogeneous and the 
radial flow is strong, <C mp^liTTQ). The data discussed 
in the present paper correspond to this latter limit. 

So if there is a significant negative value for i?^ q — 
Q in some data set, the shell of fire type of Buda-Lund 
solutions can be utilized, as a generalization of the present 
approach, as e.g. discussed also in refs. P^ll^ . Such shell 
of fire structures are also seen in blastwave models, see e.g. 
Figs. 18 and 19 of ref. |2j and also in cascade calculations 
e.g. in Fig. 12 of ref. [12]. 

The azimuthal dependence of correlation radii may be 
due to spatial deformation of the fireball as well as due to 
azimuthal dependence of the expansion velocity field. The 
interplay of these two effects in framework of the blast- 
wave model has been studied in [33] . In order to compare 
with that paper we parameterize the transverse scales with 
the help of average radius i? and asymmetry parameter 



X = asR 



Y 



R 



(40) 



(note that as corresponds to a of [33]). 

The amplitude of oscillation trivially scales with the 
absolute size of the fireball. We want to scale out this 
effect in order to directly see the effect of eccentricity. 
Thus we shall be interested in scaled oscillation amplitudes 
R^2/Rs0- (^^o^ outward radius the result would be 
similar — scaled amplitude would be just slightly smaller 
due to added f3gAT^ term in the denominator.) Such a 
ratio then does not depend on R and from eqs. ([28|) . (|36cp . 
and (j^Oj) we derive 



R 



s,2 



R'i I 



(41) 

We plot these scaled oscillation amplitudes in Figure [2] 
as functions of as and p2- The correlation between spatial 
and flow asymmetry is well pronounced: one can get the 
same oscillation with different pairs of values for the two 
parameters. Note that the dependence on flow anisotropy 
P2 in the Buda-Lund model is much stronger then in the 
blast-wave model of refs. [5] and [33 (Model 1 of the latter 
paper). 

There is an additional parameter in the Buda-Lund 
model, which influences the oscillation: the transverse tem- 
perature gradient a defined in eq. ([T2)) . We plot the de- 
pendence of scaled oscillation amplitude on a and in 
Figure [3| From eq. (|1T|) one can infer that the amplitude 
does not depend of a if p2 = 0. For the plots we have 
chosen p2 — 0.3. Increasing a tends to wash away the de- 
pendence of the amplitude on p2, as seen from eq. (j4ip . 



5 Strategy of analysis 

In this section we outline how one can infer the model pa- 
rameters from data on (i) azimuthally integrated spectra. 




-0.4 -0.2 0.2 0.4 



P2 



Fig. 2. Curves of constant Rs,2/Rs,o calculated as functions 
of as and p2 at fixed a = 0. Temperature To = 163 MeV and 
po = 0.61 (as obtained in Section [S}. Increments between the 
contours are 0.1, thick lines show -R^, 2/^3,0 — 0- PS'ir momen- 
tum Kt was chosen 0.2 GeV (upper panel) and 0.5 GeV (lower 
panel) . 



(ii) elliptic flow, (iii) azimuthally and transverse momen- 
tum dependent correlation radii. 

For a simple analysis strategy, one might observe, that 
approximate formulas can be derived assuming that the 
flow anisotropy parameter is small, p2 1- Keeping the 
leading order terms in this expansion, the following ap- 
proximative formulas are obtained: 



T,tf = n + Mpl + 0{pl), 
To 



M = mt 



V2 



To + mto^ 

P2 



1 



2m,Mp2 (1 + ^)2- 



(42) 
(43) 

(44) 



Here, Tcff is the effective slope parameter of azimuthally 
integrated mt — m spectrum . Let us observe, that this 
effective temperature depends on the freeze-out tempera- 
ture scale To and on the radial flow coefficient po-, but to 
a leading order, it does not depend on the flow anisotropy 
parameter p2. In contrast, the elliptic flow is increasing 
linearly with p2. Thus, naively, one expects that the mass 
dependence of the slope parameters can be utilized to de- 
termine the freeze-out temperature Tq and the radial flow 
parameter po, then using these parameters one could ex- 
tract the flow anisotropy parameter p2 from the data on 
the elliptic flow V2- 

If the temperature inhomogeneity parameter is very 
small, ^ To /rat, then M « mt, and one can adopt 



Mate Csanad et al.: Interplay among the azimuthally dependent HBT radii and the eUiptic flow 



7 



^<tr 



:<o: 
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PHENIX 200 GeV Au+Au spectra 20-30% (PRC69,034909 (2004)) 



0.2 0.4 0.6 0.8 1 
a 

Fig. 3. Curves of constant Rs,2/R's,o calculated as functions 
of Us and a at fixed po = 0.3. Temperature To = 163 MeV and 
po = 0.61 (as obtained in Section [5]). Increments between the 
contours are 0.1, thick lines show -Rs, 2/^3,0 — 0. Pair momen- 
tum Kt was chosen 0.2 GeV (upper panel) and 0.5 GeV (lower 
panel) . 



such a simple minded strategy. Let us then consider the 
HBT radii in this a ~ case: 



Rs 



R*.x — X 



R'*^,y — ^ 



2 

2 ^0 



+ 



r>2 _ d2 

— — cos Zip 



To + nitX^ 
To 

To + mtY^ ■ 



(45) 
(46) 

(47) 



As X and Y are given by po and p2 the new fit param- 
eters will be the actual sizes X and Y. So if a = 0, one 
can obtain the freeze-out temperature and the radial flow 
parameter from the mass dependence of the slopes, then 
momentum anisotropy parameter p2 can be obtained from 
V2 using To and po, finally X and Y come from the HBT 
radii. However, the temperature inhomogeneity parame- 
ter couples to all expressions through the effective mass 
scale M defined above, and makes such a simple strategy 
impractical, except in the case of a = 0. 

As such a naive strategy cannot be applied in the gen- 
eral case, because parameter a appears in a non-trivial 
manner in both the slope, the elliptic flow and the HBT 
data, we have tried to do a simultaneous description of 
slope parameters, elliptic flow and HBT radii. To estimate 
a reasonable range of the model parameters, we first have 
selected a centrality class, 20-30%, where such data are 
available. 



100 



% 10 




Fig. 4 

proton 



m,-m [GeV/cT 

, The exponential fits to PHENIX identified pion, kaon, 
(and antiparticle) spectra from 20-30% centrality class. 



Particle; 


m (MeV) 


Tcff (McV) 


7r+ 


139 


210.2 ± 0.8 


K+ 


494 


284.4 ± 2.2 


P 


938 


414.8 ± 7.5 


7r~ 


139 


211.9 ± 0.7 


K- 


494 


278.1 ± 2.2 


P 


938 


434.7 ± 8.6 



Table 1. Exponential slope parameters extracted from 
PHENIX weak decay corrected identifled particle spectra in 
20-30 % most central Au+Au collisions at ^/sNN = 200 GeV. 
Errors are statistical only. 



First, we have obtained the slopes of single-particle 
azimuthally integrated spectra in mj — m. The exponen- 
tial fits to PHENIX identified pion, kaon, proton (and 
antiparticle) spectra from 20-30% centrality class are ob- 
tained in just the same way as it was done for other 

,2 

centrality classes in ref. [37]: JL^..\ was fitted to 

^ ■exp(- 



{2-KratdTntdy) 

Tj-Tr — T^r — 1 — \ — J- The fit range was 0.2 GeV < 

Jtt 1 c{f{l cif+m) Jeff ' ^ 

rrit — m < 1 GeV for tt"*" and tt", while it was 0.1 GeV 
< mt — m < for , K~ , p and p. The fits are shown 
in Figure 3] and effective temperatures are summarized in 
Table [H PHENIX spectra were used here because they are 
corrected for weak decays. 

In the second step we describe qualitatively the RHIC 
20-30% central = 200GeV Au-hAu data. We com- 

pared the full Buda-Lund formulas for the slopes, ellip- 
tic flow and azimuthally sensitive HBT radii, eqs. (fT8|) . 
(HH), ((28b)) . and (pSc]) to azimuthally sensitive HBT 

radii from STAR [39], the mass dependence of effective 
temperatures from PHENIX [37|, and pion V2 data from 
STAR [3H]. For such an overall procedure, the fit parame- 
ters are To, X, Y, po, p2 and a. One can get R^.^^ and R*^y 
via eqs. ([39]) . The Buda-Lund calculations are compared 
with data in Fig. |S] for the following parameter values: 
To = 163 MeV, pn = 0.61, pa = 0.17, X^ = 33.5 fm^. 
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0.3 0.4 

p, [GeV/c] 

Fig. 5. The top panel shows the Buda-Lund calculations com- 
pared to slope parameters obtained from PHENIX weak decay 
corrected invariant momentum distribution data [37] in the 
20-30% centrality class. The middle panel shows a comparison 
to pion V2 data from STAR ^8] in the same centrality class. 
The bottom panel shows the Buda-Lund comparison to az- 
imuthally sensitive HBT radii from STAR ^39^, also measured 
in 20-30% central ^fs^ = 200GeV Au-|-Au coUisions. We used 
the full Buda-Lund formulas for the slopes, elliptic flow and 
azimuthally sensitive HBT radii, eqs. (0, (EU, (f28b)) . 

and (|28cp . and the same parameter values of To = 163 MeV, 
po = 0.61, p2 = 0.17, = 33.5 fm^, = 33.9 fm^, and 
= 0.16 were used in all the three panels. 



Y'^ — 33.9 fm^ and — 0.16. These results indicate that 
the theoretical formulas developed in the present paper 
are reasonable and could be utilized to describe qualita- 
tively the azimuthally sensitive data at in ^/snn = 200 
GeV Au-I-Au coUisions at RHIC. 



6 Conclusions 



We have shown that in the Buda-Lund model the expan- 
sion of the velocity field is oriented in such a way that the 
elliptic flow only depends on the asymmetry of the expan- 
sion velocity in the directions in and out of the reaction 
plane. Unlike the blast-wave model, there is no influence 
of the elliptic shape of transverse cross-section. 

Oscillation of the transverse correlation radii is shaped 
by flow as well as spatial anisotropy. This is again differ- 
ent from the blast-wave model: in the azimuthally non- 
symmetric generalization of that model which was able to 
describe RHIC data, correlation radii mainly depend on 
the ellipsoidal shape and less so on the anisotropy of the 
flow. 

A first analysis of the single particle spectra, elliptic 
flow and azimuthally sensitive HBT radii was presented. 
The formulas obtained from the full Buda-Lund model 
yield a reasonable description of these data with a reason- 
able choice of the model parameters, indicating, that the 
ellipsoidally symmetric Buda-Lund model has reasonable 
properties when compared to azimuthally sensitive HBT 
data. 

Note that the present formulas were obtained in a 
Gaussian approximation to the proper-time distribution 
and to the saddle-point integration of the source func- 
tion. Non-Gaussian features were reported earlier when 
performing an advanced saddle point integration - devia- 
tions from Gaussian HBT correlations were predicted in 
the longitudinal direction already in the azimuthally sym- 
metric version, see e.g. |22) . Also, a non-Gaussian proper- 
time distribution would result in non-Gaussian distribu- 
tions in directions coupled to the out direction, but these 
effects are not explored in the current manuscript. 

A unique feature, which could eventually falsify the use 
of the single saddle point approximation and the Gaus- 
sian version of the Buda-Lund model for the description 
of data is the fact that there is only explicit and no im- 
plicit azimuthal dependence of the correlation radii. In 
such models oscillation amplitudes of R^, R^, and R^^ are 
the same. Precise data could resolve if this is really the 
case. 

However, some of the conclusions like the fireball freezes 
out early, when X ~ y but X > Y, seem to be robust fea- 
tures of the data, as these features are obtained in differ- 
ent models like the Buda-Lund and the blast-wave models 
independently. 

More work has to be done to explore more quantita- 
tively the fitting procedure, the error estimation and the 
confidence level tests, using the full formulas. 
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